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We consider the phenomenology of a dimension-four operator that violates electromagnetic gauge 
invariance. Its magnitude is severely constrained by the lack of scattering of very low energy 
electromagnetic radiation off of the Cosmic Microwave Background (CMB) and by the lack of 
an induced mass when photons propagate in the CMB. We also discuss possible Lorentz-violating 
extensions of the operator basis. If a bare Proca mass exists and dominates over the induced mass, 
there is also a tight constraint from high energy scattering. 



1. INTRODUCTION 

Given that our present theories are built using the 
theme of gauge invariance, it is of some interest to under- 
stand the experimental limits on gauge invariance. More- 
over, there are motivations for possible small violations 
of gauge invariance, discussed below in Sec. 2 In this 
paper we study the effects of the Lagrangian 

C gv = -\ K A^A V A V (1) 

o 

which represents the next gauge-violating operator after 
a photon mass term. We will derive constraints of or- 
der k < 1CP 15 from the scattering of very low energy 
electromagnetic radiation, k < 10 -28 from the genera- 
tion of Lorentz-violating photon mass when propagating 
in the Cosmic Microwave Background (CMB) and yet 
stronger constraints in the possible presence of a Lorentz- 
conserving photon mass. 

The operator of Eq. 1 leads to a cross section for 
photon-photon scattering that diverges at low energy. 
This can be seen simply from the dimensionless coupling 
constant k and the absence of a mass such that the total 
cross section must go as 



simply on dimensional grounds. Here s = E 2 cm is the 
square of the center of mass energy. Such an inter- 
action would scatter the lowest energy electromagnetic 
radiation. In particular, radio waves from distant ob- 
jects would scatter off the Cosmic Microwave Background 
(CMB) and would directly not reach us if the mean free 
path is too short. We study this effect in Sec 3 and use it 
to bound k. Section 4 considers some Lorentz violating 
variants of this interaction. 

Moreover, when propagating through the Cosmic Mi- 
crowave Background (CMB) a photon will pick up an ef- 
fective mass term from this interaction. This mass is not 
Lorentz invariant as it depends on the rest frame of the 
CMB. Somewhat remarkably, such a Lorentz-violating 
mass is better behaved than a Lorentz-invariant one, 
and has interesting phenomenology which we explore in 



Sec. 5. Finally we explore the possible enhancement of 
the limits that happens for some values of a Lorentz- 
conserving mass term, which occur because scattering 
from transverse to longitudinal gauge bosons occurs with 
a cross section inversely proportional to the possible pho- 
ton mass. 



2. BRIEF COMMENTS ON GAUGE 
NON-INVARIANCE 

In looking for signals of new physics, studies of vio- 
lations of symmetries are particularly useful. Much of 
the focus of searches for physics beyond the standard 
model uses effective Lagrangians of dimension greater 
than four, which preserve the gauge symmetries of the 
Standard Model, but allow the violation of any global 
symmetry Another significant sub-field looks for vi- 
olations of Lorentz invariance. This involves a further 
extension of the Standard Model [2], this time involv- 
ing Lorentz violating operators of dimension 2, 3, 4 etc, 
still with the assumption that the Standard Model gauge 
symmetries are exact. Here we are concerned with a piece 
of phenomenology giving up the last of the symmetries - 
gauge symmetry. 

One motivation comes from the connection between 
Lorentz invariance and gauge invariance. It is a little ap- 
preciated feature that the massless gauge fields are not 
themselves Lorentz covariant|3j. The Lorentz transfor- 
mation needs to be supplemented by a gauge transfor- 
mation in order to make an overall invariance. A simple 
example is a Lorentz boost of a transverse photon polar- 
ization vector e^(k) in a direction that is not along the 
initial propagation direction k v , resulting in a polariza- 
tion vector which is no longer transverse to the propaga- 
tion direction in the new frame. A gauge transformation 
is required to make this polarization vector transverse. 
Because of this connection it is possible that loop ef- 
fects with Lorentz violating interactions could produce 
gauge non- invariant effects. Indeed such effects are found 
in specific calculations There are some controver- 
sies concerning the regularizat ion of the gauge- violating 
effects ||j but a reasonable conclusion could be that the 
generation of such effects is possible depending on the 



2 



nature of the ultimate high energy theory. 

Another motivation is the possibility of emergent 
symmetries There are known examples in condensed 
matter systems where a gauge symmetry emerges in the 
ground state of a theory that does not originally have 
such a symmetry. There are also many attempts at emer- 
gent space-time, which would require that the general co- 
variance of gravity be emergent also. In such cases, the 
leading approximation to the full theory would satisfy the 
emergent gauge symmetry. However, since the underly- 
ing theory does not respect the symmetry, there can be 
small corrections to the leading approximation violating 
the symmetry. At the very least, loop diagrams should 
be able to probe the lack of symmetry at the fundamental 
scale, as the loops involve integrations over all energies. 
Looking for violations of gauge invariance is the part of 
the phenomenology of emergent gaug e sy mmetry 11 1 . 

The Weinberg- Witten theorem [ljj LL3| connects these 
two motivations as it says that it is not possible to have 
a composite or emergent gauge boson 1 or graviton from 
a Lorentz invariant theory. This implies that emergent 
gauge symmetry would come from a Lorentz violating 
theory, and Lorentz symmetry itself could be emergent. 

The theorem does not give any indication of how large 
any residual violations are, but it make sense to search 
for gauge symmetry violation and Lorentz violation to- 
gether. Because of this connection we will also briefly 
study a generalization of Eq. 1, 

C = -\f 2 - X -K^ va ^A v A a Af i (3) 

with a more general tensor K^ap. Eq. 1 is reproduced 
with 



(4) 



and we will consider the case where there is a preferred 
direction in space b^, with the possibilities 



and 



KpvaP — n'b^bvbabp 



(5) 



(6) 



only at the highest energies, leading the the well-known 
GZK cutoff of ultra- high-energy cosmic raysjllj]. How- 
ever, with the interaction of Eq. 1, the strongest scatter- 
ing is at low energy. 

Let us calculate the cross-section for 77 — > 77 scatter- 
ing using Eq. 1. The matrix element is found to be 

M = fc[e(pi,Ai) • e(p 2 , A 2 )e(p 3 , A 3 ) • e(p 4 ,A 4 ) 
+e(pi, Ai) • e(p 3 , A 3 )e(p 2 , A 2 ) ■ e(p 4 , A 4 ) 
+e(pi, Ai) • e(p 4 , A 4 )e(p 2 , A 2 ) • e(p 3 , A3)] (7) 

At this stage we do not consider a photon mass term, 
returning to this topic in later sections. Without a mass, 
we consider only the production and scattering of trans- 
verse photons. Operationally, this means the polarization 
vectors are 'physical', 



e°(p,±l)=0 
p- e{p,±l) = 



(8) 
(9) 



Because the interaction is not gauge invariant, the usual 
replacement Y, polarizaUon e M e v -> -g^ can not be made 
here. This simple replacement only works in QED be- 
cause of gauge-invariance leading to the cancellation of 
the contributions coming from the unphysical longitudi- 
nal and scalar degrees of freedom. To square the invari- 
ant amplitude then, we directly perform the calculation 
in the CM frame for each polarization configuration in- 
dividually. To this end, we employ a linear polarization 
basis with the following convention, 

e(p,l)xe(p,2)=p (10) 

The non-vanishing transitions are shown in Table I. We 

TABLE I: Polarizations and Matrix Elements 



Config. 



M 



12 
22 
11 
22 
21 
21 
12 
12 



11 
22 
22 
11 
21 
12 
21 



K COS 9cm 
K COS 9 cm 
K COS 9cm 
K COS 9cm 

3 k 

K 
K 



12 K(l + 2cos 2 6> em ) 



3. SCATTERING OFF THE COSMIC 
MICROWAVE BACKGROUND 

Photons from a distant source can scatter off of the 
CMB. With the usual gauge invariant interaction, the 
Euler Heisenberg Lagrangian, this becomes significant 



1 The theorem applies to Yang-Mills gauge bosons but not to those 
of a U(l) symmetry. However, in the Standard Model the phys- 
ical photon is a mixture of SU(2) and U(l) fields so that we 
expect that the theorem is applicable to real photons. 



sum (average) over the initial (final) polarization states 
and obtain the modulus-squared amplitude, 

i \M\ 2 = n 2 {3 + 2cos 2 8 cm + cos 4 9 cm ) (11) 

polarizaton 

where 9 cm is the angle the outgoing particles makes with 
the collision axis. The differential scattering cross-section 
is readily found to be, 



da 



dVt cm 6An 2 s 



(3 + 2 COS 2 9cm + cos 4 9 C 



(12) 
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which upon integration yields the total cross-section, 



29k 2 

1207TS 



(13) 



In order to obtain a bound on the coupling constant n, 
we consider the scattering of a low energy photon (E 7 ) 
cmmittcd from a distant quasar off a CMB photon (E) 
and calculate the collision rate for such a reaction. Obvi- 
ously we perform the calculation in the earth frame since 
all energy measurements are known only with respect to 
the aforementioned frame. We align the momentum of 
the incoming photon along the z-axis and consider the 
momentum of the CMB photon to make an angle 9 with 
the z-axis. Thus we have, 



s = 2EE 1 (l - cosf 
= 4EE~, sin 2 - 



(14) 



The CMB photons can impinge on the incoming photon 
at any angle and any energy. Moreover, they obey Bosc- 
Einstcin statistics and we average appropriately over all 
angles and energies. The total cross-section diverges be- 
cause of the pole at 9 = 0° . This comes from the limit 
of collinear photons. However, we will see that the total 
collision rate is appropriately finite. 

The total collision rate, whose inverse yields the mean 
free path, can be used to bound k. Following [la ], we 
define the number density of CMB photons as 



dp(E,9,< 



E 2 dEdcf)dO sin 6 
4^3 (exp/3E - 1) 



(15) 



The rate of collision between the incoming photon and 
CMB photons becomes 

dT = \v~, ~v E \adp(E, 6,<f) (16) 

where |u 7 — ve\ is the relative speed in the earth frame. 

v 1 — ve = — (sin 9 cos </>, sin # sine/), cos 9—1) (17) 



ItL — Ve\ =2 sin - 



(18) 



Finally we integrate over all angle and energies to obtain 
the total collision rate, 



oq K 2 

T = ? / dEdcos9 



E 2 sin 



12(),T :; ./ 4E£ 7 sin 2 § (ef> E -l) 

29n 2 ((2) 



r = 



Here 



120tt 3 £L/3 2 



C(2) = 1.645 



(19) 
(20) 



(21) 
(22) 



It is well known that very low-energy radio signals 
reach the earth coming from distant quasars, which are 



billions of light-years away from us. We place an upper 
bound on n by demanding that the time taken by these 
low-energy signals to reach us has to be less than the 
time between individual collisions, given by the inverse 
of the collision rate. 

The most distance known radio-loud quasar is the 
source J1427+3312 [Hj]. It has been detected at 1.6 GHz 
with a redshift of z=6.12. The light-travel distance is 
simply calculated from the redshift, and thus the light- 
travel time is 



t « 4 x 10 17 sec 



(23) 



We also work with the average, 



</? 2 > = 



(2.725) 2 (z + 1) 



2.547 x 10 6 (eVT 2 (24) 



In the above, we used Tcmb — 2.725(1 + z). This leads 
to the bound 



k < 3.686 x 10 



-15 



(25) 



In fact we can do slightly better by using a closer source 
which is however observed at a lower frequency. For ex- 
ample in reference [ljj , the sources in the sample had ob- 
servations in the frequency as low as 74 MHz. The high- 
est red-shift among the radio-loud quasars is the source 
3C 446 with z=1.404. We can use the data of this source 
to find 



T « 2.857 x Wsec 



Then the final bound is 



{P 2 



7.556 x 10 6 (eVT 2 

(26) 



k < 1.615 x 10 



-15 



(27) 



We see that the trade-off between energy and red-shift 
produces an 0(1) effect. Finally, we note that it is very 
likely that space-based observations could detect quasars 
in lower regimes of the radio spectrum which greatly en- 
hances our bound. 



4. COMBINED LORENTZ AND GAUGE 
VIOLATION 

Because of the possible connection of Lorentz and 
gauge violation, we here consider the combined violation 
of both types of symmetries. The resulting bounds are 
not much different, but there are some interesting fea- 
tures which emerge when both symmetries are broken. 

Lorentz violation is generally parameterized by a fixed 
background vector (or tensor) that specifies a preferred 
direction [2j. For our studies we use a particular preferred 
frame designated by b^. In this case the preferred frame 
can also enter the gauge violating interaction. 

We first consider the following coupling, 



^ /if a/3 



= n'b^bvbabp 



(28) 
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If we again consider the scattering of transverse photons 
with this modified coupling, we find for the matrix ele- 
ment 

M = 3k' [b- e(pi, Ai) b-e(p 2 ,A 2 ) b-e(p 3 ,X 3 ) b ■ e(p 4 , 

(29) 

To compute the cross section, we carry out the polar- 
ization sum in the CM frame. To this end, we use the 
following polarization sum formula Q 



The squared matrix element is then, 

M 2 = ^' 2 (P-bl) 2 ^-(b-p 3 y 



(30) 



(31) 



Where 63 is the component of b along the collision axis in 
the CM frame. Notice that the differential cross section 
is not azimuthally symmetric in this case as the back- 
ground vector field defines a preferred direction in any 
single frame. The cross section is then, 



is not invariant under the so-called 'observer' Lorentz 
transformations which is just the conventional class 
of transformations of special relativity. As discussed 
before, the polarization vectors of the gauge-field are 
A 4)kot proper four vectors, thus the scalar products in the 
matrix element are not strictly Lorentz invariant. This 
is different than the framework of the Lorentz violating 
SME, which is invariant under the aforementioned class 
of transformations but not invariant under the so called 
'particle' Lorentz transformations, where breaking the 
symmetry comes about by the existence of background 
constant tensor fields which define a preferred direction 
in any single frame. As shown, adding gauge-violation 
to the framework of SME adds this new feature which 
enables us to probe Lorentz-symmetry violation in the 
conventional sense. A final remark worth mentioning is 
th at the pure gauge-violating theory does not exhibit 
such a feature because the matrix element includes 
only scalar products of polarization vectors. In other 
words, the field is purely self-interacting so that Lorentz- 
invariance violation is 'hidden' but in the present case, 
the background vector acts like an external field coupled 
to the gauge-field and thus the result we have just shown. 



80ns 



(32) 



Let us for completeness also consider the following cou- 
pling, 



(33) 



As before, k" is the parameter controlling the strength 
of gauge-invariance violation, while 6 M is a constant four- 
vector field controlling Lorentz-invariance violation. 

We follow the same procedure as last section and com- 
pute the matrix element for each polarization configura- 
tion individually in the CM frame. As the matrix ele- 
ment for all polarization configurations are non-trivial, 
we choose to compute the cross section for a specific con- 
figuration, namely .M(ll — > 11) where we again used a 
linear polarization basis with the conventions of the last 
section. 



cr(ll 11) 



„."2 



15360tts 
+ 142& 2 fr 



(5096? + 
432fr 2 &; 



9b A 2 



- 646 4 

mlbl) 



(34) 



The inclusion of Lorentz violation does not generate 
any improvement in the bounds. The constraints from 
the mean free path go through as before. We neglect 
the potential dependence on the relative direction of a 
particular quasar to the spatial direction of b and then, 
taking into account the differing numerical factors in the 
cross sections, quote rough bounds 



k'6 4 < 10 



-14 



n"b 2 < 1CT 14 



(35) 



We note that the combined Lorentz and gauge viola- 
tion posesses a very interesting feature, namely that it 



5. MASS GENERATION FROM THE CMB 

Once one considers a gauge non-invariant interaction 
such as Eq. 1, there is no symmetry forbidding a photon 
mass. However, the interaction itself does not have to 
generate a photon mass in a vacuum. A loop diagram in- 
volving two of the photons of the 4 photon vertex would 
have the potential to generate a mass term. However, 
the loop integral vanishes when regularized dimension- 
ally. Despite this, the interaction of Eq. 1 will necessar- 
ily generate a mass for a photon moving in a background 
of other photons, i.e. the CMB. We evaluate this mass in 
this section and explore some of its more unusual prop- 
erties. 

We treat the heatbath as a background field which the 
photon is propagating in. This method is equivalent to a 
photon loop in real-time finite temperature field theory. 
We treat the interaction Lagrangian by splitting the field 
into a background piece - the heat bath - and a quantum 
propagating photon, 

(36) 



A^ 1 = Al + A^q 



Expanding the product of fields in the Lagrangian, we 
collect the six terms quadratic in the background field. 

nquad _ K ^aP ( All AV Ad A0 i A^ A a A v A$ 
— o \ b b Q Q ' A b b Q O. 



*b ^-b ^Q^Q ^ ^ h ^ b ^Q^Q 
j- A^ A@ A a A" 4- A a A v A^ A@ 
-lAP A v A a A^ 4- A a A 13 A^ A v 



H ^b ^Q^Q) 

The finite-tempreture ground state is defined as, 
(/3\ a i(p,\)a(q,\'M =n B (E ff )(2ir) 3 S X y5 3 (p-q) (37) 



■5 



Where ns is the usual Bosc-Einstcin distribution for pho- 
tons. The contribution of the CMB would be calculated 
by the taking the ground state expectation value of the 
quadratic Lagrangian. First, we evaluate 

{P\AW\fi)=Y, J j^- p n B (E p )[e^p,X)e v (p,X)} 

(38) 

To evaluate the above integral, we note that from symme- 
try the tensor integral is proportional to the Kronecker 
delta. 

£ / j0E; nB{Ep) t £i(p ' A)e ' (p ' A) ] = ^ (39) 



Contracting both sides with Sij, we find 



Where I is a simple integral over the BE distribution, 



Thus the quadratic Lagrangian reads, 
1 / 2 

£quad _ r / /1 2 



2 Kl { Q ~ 3 A Q' A Q 



(42) 



Now we construct the free field Lagrangian for the quan- 
tum field, which reads 



£ = -^F 2 + ±m a pA a AP 



Where we defined, 



moo = k! rriij = 



(43) 



(44) 



We see that k has to be positive definite to ensure the 
theory does not have growing exponential solutions. 
In order to solve the equations of motion, we consider 



(45) 
(46) 



By acting d v on the equation of motion and requiring 
the external current coupled to the field to be conserved, 
we obtain the following constraint which ensures the field 
has only three degrees of freedom. 



d v d ll Fi»' + m v a (d v A a ) = Q 
m v a (d v A a ) = 

d A° + ^A* = 
We propose the usual wave ansatz, 

A(x) oc e -<(<"*-**) 



(47) 
(48) 

(49) 



(50) 



Thus Eq. 49 becomes, 
3 



luA° - k ■ A = 



(51) 



Due to the manifest Lorentz non-invaraince, the quanti- 
zation procedure needs some care. First, we explore the 
dispersion relations by plugging the wave ansatz in the 
equations of motion. 

(-L0 2 + fc 2 ) A v + k u (k- A)+m v a A a = (52) 

For definiteness, we take k = kz. Taking v — and using 
the constraint equation, we obtain 



(-u? + fc 2 ) A + (uA° - 



» + m Q A° = (53) 



(40) This yields the following dispersion relation 



5" 



2 k 2 = Kl 



(54) 



The same dispersion relation is readily obtained if v = 3. 
On the other hand, for v = 1, 2 we get 



2 + &)A 1 ' 2 + UlA 1 * = 



which yields the dispersion relation, 



oj 2 - k 2 = Ul 



(55) 



(56) 



The differing dispersion relations are a clear manifesta- 
tion of the Lorentz non-invariance of the theory. We next 
show how to quantize the theory by introducing differ- 
ent dispersion relations for transverse and longitudinal 
modes. 

For later convenience, we define 



-kI = m 2 



(57) 



The aforementioned theory describes three degrees of 
freedom. Guided by the previous analysis , we proceed 
by proposing two different wave ansatz for the different 
modes of propagation. 

2 

A k (x) oc e'^-^eiX, k) + e"'^'-^) e(3, k) 

(58) 



A=l 



where four- vectors are displayed without arrows. To sat- 
isfy the equations of motion and the constraint Eq. 51, 
the polarization vectors have to satisfy 

e°(l,fc) = e°(2,fc) = (59) 
e{l,k)-k = e{2,k)-k = (60) 

e " (3 '* )= ( 1 '^*) (61) 
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The third polarization vector shows a very interesting 
feature; it becomes null-like at some value of the 3- 
momentum and hence its normalization cannot be fixed. 
This is another manifestation of breaking Lorentz sym- 
metry. 

The two dispersion relations are, 



(62) 
(63) 



The fact that Eq. 62 has the canonical form of energy- 
momentum relation and that it is associated with the 
transerve modes is tempting to interpret m as the 'phys- 
ical' mass. 

We can employ this generated mass to set an upper 
bound on k by using the most reliable experimental upper 
bound on the photon mass. We use the CMB current 
temperature to compute /, 



I(T CMB = 2.71) = 45 x 10- w (eV) 
The photon mass is bounded to be (Til] 
m 2 ph <lx 10- 36 (eV) 2 
leading to a bound on k 



k < 1.3 x 10 



-28 



(64) 



(65) 



(66) 



6. MASS EFFECTS IN SCATTERING 
AMPLITUDES 

6.1. Lorentz invariant mass 

In this section we explore a possible scenario which 
greatly enhances our bound. The logic is that if gauge- 
violation is admissible we do not have any theoretical ar- 
gument upon which to set the mass of the photon to zero. 
Rather we should use the massive Proca theory, treat the 
photon mass as a parameter and assign its value from ex- 
periment. In this case, there is a longitudinal polarization 
and we can again consider scattering from transverse to 
longitudinal polarizations. In this situation however, the 
cross-section goes as to -4 and becomes extremely strong 
for low masses. For many values of the mass, this leads 
to a very tight constraint on n. 

The Lagrangian is, 



C = — F + — A - -A 4 
4 2 8 



(67) 



Much like Section 3, we look at the process where now 
massive photons traveling extra-galactic distances scat- 
ter off CMB photons. A few remarks are worth men- 
tioning at this stage regarding our choice of considering 
transverse photons in the initial state. Massive vector 



bosons are spin-1 particles having three degrees of free- 
dom represented by polarization four-vectors satisfying 
the following relations 



e(k,X)-k = , e(fc,A)-e(fc,A') = -*AA' 
^2 e M (fc, A)e„(fc, A) = -g^ 



(68) 
(69) 



spin 



In any frame of reference other than the rest frame, the 
polarization vectors can be conveniently chosen as fol- 
lows; two are taken to be purely spatial and transverse to 
the 3-momentum following a right-hand rule. If we work 
with states corresponding to circular polarization, then 
the transverse vectors correspond to the z-projection of 
the spin being ±1. 



e(fc,l) ■ k = e(k,2) ■ k = 
e(fc,l) x e{k,2) = ft 



(70) 
(71) 



While the third vector is longitudinal to the 3- 
momentum, which corresponds to the z-projection of the 
spin being 0. 



e"(ife,3) = — (Iftl.wft 
m 



(72) 



In a frame where u> 3> m, we can expand the components 
of the longitudinal vector and obtain 



e"(fc,3) 



ft' 1 



O 



(73) 



We imagine those massive photons being produced at 
extra-galactic objects via coupling to a conserved elec- 
tromagnetic cuurent and thus, the Ward identity holds. 
Therefore, given the expansion of the longitudinal polar- 
ization vector at high energies, any matrix element with 
an on-shell longitudinal photon vanishes to leading or- 
der in the energy. In an arbitrary frame where u) S> m, 
the polarization of a photon with four-momentum ft is 
generally a superposition of the transverse and longitu- 
dinal polarizations, but based on the last discussion the 
contribution from the longitudinal mode vanishes. 

Following the previous comments, we calculate the 
cross-section for the scattering of two massive transver- 
sally polarized photons in the earth frame. For simplicity, 
we take the initial state to be a head-on collision aligned 
on the z-axis and, without loss of generality, we work 
with linear polarization basis. Accordingly, the initial 
state polarization vectors are 



e"(pi,l) = (0,1,0,0) 
c*(fc,l) = (0,0,1,0) 



(74) 



The fact that the collision is head-on simplifies the com- 
putation considerably, the CM frame is related to the 
earth frame via a simple boost along the z-axis with the 
boost factor found to be, 



u>\ - uj 2 



(75) 
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This boost obviously does not affect the polarization vec- 
tors of the incoming photons and moreover we average 
over the final state polarizations using the spin sum for- 
mula Eq. 69. The squared matrix element is found to 
be, 



9 ^ 



9m 4 



[(m 2 + (p 3 ) 2 )(m 2 + (p 2 ) 2 ) (76) 



+ (m 2 + { P IY) (m 2 + (p\f) + 2 P l P l P \ P l] 

All momenta are measured in the CM frame. 
We consider the high-energy limit of the reaction, and 
therefore we neglect the masses compared to the mo- 
menta. The high energy behavior is dominated by the 
production of longitudinal photons in the final state. The 
cross-section is then found to be, 



1 



d 3 p 3 d 3 p 4 



4|wi|s J (2tt) 3 2w 3 (27r) 3 2w 4 
{2n) A 8{p 1 +p 2 -p 3 -p i )\M\ 2 

We carry out the integral over d 3 pi to find, 



(77) 



11527r 2 |wi|m 4 s 



d 3 P3 



*(w3-Wi)((pS) 2 + (p1) s 



(78) 



Performing the phase space integral and expressing the 
result in terms of the Mandelstam variable, we get 



(s - 4m 2 ) 2 



21607TTO 4 



(79) 



We note that if we had taken longitudinal photons in the 
initial state, the resulting cross section would go as m~ 8 . 
However, we do not use this in producing constraints, 
because we do not expect to have longitudinal photons 
produced in quasar emission, and would not miss them 
if they scattered before reaching us. 

We place a bound on n by demanding the mean free 
path between collisions is greater than the distance trav- 
eled by the photon from the source to earth. Looking at 
the cross-section, the tightest bound is obtained by con- 
sidering the scattering of very energetic photons, namely 
GRB signal. We use 10 5 eV, which is the typical mini- 
mum energy in GRB's, and work with the mean CMB 
energy ojcmb = 6.34 x 10~ 4 eV7 The mean free path is 
given by[14| 



A = (n g a) 



(80) 



Here, n g is the mean number density of CMB photons. 
Although we only considered head-on collisions, we use 
n g = 410.4cm -3 . The effect of this approximation is ob- 
viously of 0(1) as the GRB energy could be well above 
\Q 5 eV . The typical distance of GRB sources is billions 
of light years, L ~ 10 27 cm. Finally, we use the upper ex- 
perimental bound on the Lorentz-invariant photon mass 
and find 



K < 0.67 x 10~ 46 



(81) 



As discussed in the next subsection and in the summary, 
this bound is only realized if the Proca mass is dominant 
over the Lorentz violating mass that comes from the in- 
teractions with the CMB. 



6.2. Lorentz- Violating Mass 

We now also compute the cross section of the same 
process but with the Lorentz-violating massive theory de- 
scribed in Sec 5. We have just seen that if the Proca mass 
is dominant, there is a bad high energy behavior from the 
scattering into the longitudinal polarization state. This 
comes from the m _1 factor seen in Eqs. l72l and l73l For a 
Lorentz violating mass, this feature does not occur, with 
the inverse mass being replaced by an inverse momentum. 

The initial state polarization vectors are, 



\P2, 1) 



(0,1,0,0) 
(0,0,1,0) 



(82) 
(83) 



The computation is carried in the CMB frame since this is 
the preferred frame in which the theory is defined. To this 
end, the final state photons are longitudinally polarized, 



e"(ft,3) = 1, 



-Pa 



(84) 
(85) 



If the initial state is a head-on collision, then the matrix 
element is found to be 



M 



9KU)L,3WL,4 1 i 2 
{P3P4 



pIpI) 



(86) 



To simplify the evaluation of the phase space integral, 
we proceed by taking Wx = 0J2 exploiting the freedom to 
choose the energy of the interacting photons, since the 
universe is filled with photons propagating at arbitrary 
energies over all the electromagnetic spectrum. 



The squared matrix element is then, 



, 324k 2 w£ 3 , 
= 625; p4 sin 3 sin 2 03 cos 2 <f> 3 



(87) 



Peforming the phase space integral to find the total cross 
section, 

9\/3k 2 s 



4000tt ( s _ 4 TO 2)^ 

In contrast to the Proca case, this cross section is well 
behaved at high energy. The Lorentz violating nature of 
the third polarization state softens the high energy be- 
havior. We therefore do not obtain any extra constraint 
from high energy scattering if the Lorentz violating mass 
is dominant over the Proca mass. 
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FIG. 1: Domain of k and the Proca mass m p . As explained 
in the text, regions I and II are allowed and regions III, IV, 
and V are forbidden 

7. SUMMARY OF BOUNDS 

Figure [1] shows the allowed domain of k, which is en- 
tangled with the relative size of the possible Proca mass. 
Let us explain the constraints given here. Regions IV and 
V are excluded because the photon mass is too large. For 
region V, this is the standard exclusion region of a Proca 
mass. Region IV corresponds to the constraint from gen- 
erating a Lorentz violating mass through the interactions 



with the CMB, as described in Sec 5. For combinations of 
k and the Proca mass that satisfy the above constraint, 
the situation depends on the relative size of the Proca 
mass versus the Lorentz violating mass from interactions 
with the CMB. If the Proca mass is dominant, there is 
the constraint from high energy scattering into the longi- 
tudinal polarization state, described in Sec 6.1, that rules 
out region III. However, if the Lorentz violating mass is 
dominant, then the longitudinal polarization state is well 
behaved at high energy and there is no such constraint. 
This yields region I as an allowed region. Finally in re- 
gion II, k is so small that it satisfies all constraints. 

The interactions of the photons thus provide extremely 
strong constraints on the dimensionless parameter k that 
governs the strength of this gauge violating interaction. 
This of course is consistent with the standard expecta- 
tion that such an operator is forbidden by gauge invari- 
ance. Combined with constraints on the photon mass, 
our work helps quantify the degree that this expecta- 
tion is supported by evidence. Of course, it is possible 
that gauge violation might only show up in higher di- 
mensional operators. The mass term carries dimension 
two and our operator is dimension four, both of which 
are the dimensions allowed for renormalization interac- 
tions. If the gauge symmetry is emergent, it is plausible 
that at dimension two and four one obtains exactly the 
gauge invariant theory and that the signal for gauge vi- 
olation will be higher dimensional operators suppressed 
by powers of some heavy scale. Such operators are also 
worthy of study, although the constraints are expected 
to be significantly weaker. 
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